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Abstract 

The Tutte polynomial of a graph G is a two-variable polynomial 
T(G; x, y) that encodes many interesting properties of the graph. We 
study the complexity of the following problem, for rationals x and y: 
take as input a graph G, and output a value which is a good approx- 
imation to T(G;x,y). Jaeger, Vertigan and Welsh have completely 
mapped the complexity of exactly computing the Tutte polynomial. 
They have shown that this is #P-hard, except along the hyperbola 
(x — l)(y — 1) = 1 and at four special points. We are interested in de- 
termining for which points (x, y) there is a fully polynomial randomised 
approximation scheme (FPRAS) for T(G; x,y). Under the assumption 
RP 7^ NP, we prove that there is no FPRAS at (x,y) if (x,y) is is 
in one of the half-planes x < — 1 or y < - 1 (excluding the easy-to- 
compute cases mentioned above). Two exceptions to this result are 
the half- line x < — 1, y = 1 (which is still open) and the portion of the 
hyperbola (x — l)(y — 1) = 2 corresponding to y < —1 which we show 
to be equivalent in difficulty to approximately counting perfect match- 
ings. We give further intractability results for (x,y) in the vicinity of 
the origin. A corollary of our results is that, under the assumption 
RP ^ NP, there is no FPRAS at the point (x,y) = (0, 1 - A) when 
A > 2 is a positive integer. Thus, there is no FPRAS for counting 
nowhere-zero A flows for A > 2. This is an interesting consequence of 
our work since the corresponding decision problem is in P for example 
for A = 6. Although our main concern is to distinguish regions of the 
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Tutte plane that admit an FPRAS from those that do not, we also 
note that the latter regions exhibit different levels of intractability. At 
certain points (x,y), for example the integer points on the a;-axis, or 
any point in the positive quadrant, there is a randomised approxima- 
tion scheme for T(G; x, y) that runs in polynomial time using an oracle 
for an NP predicate. On the other hand, we identify a region of points 
(x, y) at which even approximating T(G; x, y) is as hard as #P. 

1 Summary of results 

The Tutte polynomial of a graph G (see Section [2TT]) is a two- variable poly- 
nomial T(G; x, y) that encodes many interesting properties of the graph. We 
mention only some of these properties here, as a longer and more detailed 
list can be found in Welsh's book [2*2*] . 

• T(G; 1, 1) counts the number of spanning trees of a connected graph 
G. 

• T(G; 2, 1) counts the number of forests in G (the number of edge sub- 
sets that contain no cycles). 

• T(G; 1,2) counts the number of edge subsets that are connected and 
span G. 

• T(G; 2, 0) counts the number of acyclic orientations of G. 

• The chromatic polynomial P(G; A) of a graph G with n vertices, m 
edges and k connected components is given by 

P(G; A) = (-l) n ^A fe T(G; 1 - A, 0). 

When A is a positive integer, P(G; A) counts the proper A-colourings 
of G. 

• The flow polynomial F(G; A) is given by 

F(G; A) = (-l) m - n+k T(G; 0, 1 - A). 

When A is a positive integer, F(G; A) counts the nowhere-zero A-flows 
of G. 

• The all-terminal reliability polynomial R(G;p) is given by 

R(G;p) = (l-p) m - n+k p n - k T(G;l,l/(l-p)). 

When G is connected and each edge is independently "open" with 
probability p, R(G,p) is the probability that there is a path between 
every pair of vertices of G. 
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• For every positive integer q, the Tutte polynomial along the hyperbola 
H q given by {x — l)(y — 1) = q corresponds to the partition function 
of the (/-state Potts model. 

We study the complexity of the following problem, for rationals x and y: 
take as input a graph G, and output a value which is a good approximation to 
T(G; x, y). Jaeger, Vertigan and Welsh [TU] (see Section |2~4"[) have completely 
mapped the complexity of exactly computing the Tutte polynomial. They 
have shown that this is #P-hard, except along the hyperbola H\ and at the 
four special points (x,y) G {(1, 1), (0, -1), (-1, 0), (-1, -1)}. (#P is the 
analogue, for counting problems, of the more familiar class NP of decision 
problems.) 

We are interested in determining for which points (x, y) there is a fully 
polynomial randomised approximation scheme (FPRAS) for T(G;x,y). An 
FPRAS is a polynomial-time randomised approximation algorithm achieving 
arbitrarily small relative error. Precise definitions of FPRAS, #P, and other 
complexity-theoretic terminology will be provided in Section 12.21 

It is known that there is an FPRAS for every point (x, y) on the hy- 
perbola H2 with y > 1 — this follows from the Ising result of Jerrum and 
Sinclair [12]. No other general FPRAS results are known. A few negative 
results are known — see Section 12.41 

Our goal is to map the Tutte plane in terms of FPRASability as com- 
pletely as possible. The specific contribution of this article is a substantial 
widening of the region known to be non-FPRASable. 

Our contributions are summarised in Figured) In particular, under the 
assumption RP 7^ NP, we prove the following. 

(1) If x < — 1 and (x,y) is not on Hq or Hi, then there is no FPRAS at 
(x,y) (Corollary 

(2) If y < — 1 and (x,y) is not on Hi or H2, then there is no FPRAS at 
(x,y) (Corollary [5] when (x,y) is not on Hq and Lemma [6] for the case 
in which (x,y) is on Hq). 

(3) If (x, y) is on H2 and y < — 1 then approximating T{G; x, y) is equivalent 
in difficulty to approximately counting perfect matchings (Lemma [7J). 

(4) If (x, y) is not on Hi and is in the vicinity of the origin in the sense that 
I a; I < 1 and \y\ < 1 and is in the triangle given by y < — 1 — 2x then 
there is no FPRAS (Lemma Ej). 

(5) If (x,y) is not on Hi and is in the vicinity of the origin and is in the 
triangle given by x < — 1 — 2y then there is no FPRAS (Lemma [9]). 

(6) The two previous intractability results (results (4) and (5)) can be par- 
tially extended to the boundary of the triangles (Lemma [101 and [TU . 
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(7) If (x,y) is in the vicinity of the origin and q = {x — l)(y — 1) > 1.5 
then there is no FPRAS (excluding the special points at which exact 
computation is possible) (Lemma [T2|) . 




Figure 1: Green points are FPRASable, red points are equivalent to perfect 
matchings and gray points are not FPRASable unless RP=NP. We don't 
know about white points. The points depicted in black are at least as hard 
as gray and are presumably harder — this is the region of q = 4 with 
y £ (—1,0) and approximating Tutte is actually #P-hard here. (There are 
presumably more such points.) 

Result (2) above implies that, under the assumption RP ^ NP, there is 
no FPRAS at the point (x, y) = (0, 1 — A) when A > 2 is a positive integer. 
Thus, there is no FPRAS for counting nowhere-zero A flows for A > 2. This 
is an interesting consequence of our work since Seymour [16] has shown that 
the corresponding decision problem is in P for A = 6. In particular, a graph 
has a 6-flow if and only if it has no bridge (cut edge). 

Although our main concern is to distinguish regions of the Tutte plane 
that admit an FPRAS from those that do not, we also note that the latter 
regions exhibit different levels of intractability. At certain points (x,y), for 
example the integer points on the x-axis, or any point in the positive quad- 
rant, there is a randomised approximation scheme for Tutte(x, y) that runs 
in polynomial time using an oracle for an NP predicate. On the other hand, 
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Theorem [T7] identifies a region of points (x,y) at which even approximating 
Tutte(x, y) is as hard as #P. These two kinds of intractability are very 
different, assuming #P is a "much bigger" class than NP. 

2 Definitions and context 

2.1 The Tutte polynomial 

The Tutte polynomial of a graph G = (V, E) is the two- variable polynomial 

T(G; X ,y)=Y J (x- l)^)-^)(y - !) W-n+KW (1) 
ACE 

where k(.A) denotes the number of connected components of the graph (V, A) 
and n = |V|. Following the usual convention for the Tutte polynomial [T7] 
a graph is allowed to have loops and/or multiple edges, and we use the term 
"graph" in this way except where we explicitly state otherwise. The Tutte 
polynomial is sometimes referred to as the "Whitney- Tutte" polynomial, or 
the "dichromatic polynomial". See [20] and |22j . 

2.2 The complexity of counting and approximate counting 

We start with a brief summary of the complexity of counting. See for 
more details. A counting problem can be viewed as a function / : T,* — > N 
mapping an encoding of a problem instance (encoded as a word in a finite 
alphabet, X) to a natural number. For example, / might map an encoding 
of a graph G to the number of independent sets of G. #P is the analogue 
of NP for counting problems. A counting problem / : X* — * N is in #P if 
there is a polynomial-time predicate x '■ ^* x ^* ~> {0, 1} and a polynomial 
p such that, for all instances x G X*, 

/(x) = e X* | x(aJ 3 w) A \w\ < p{\x\)}\. 

It is straightforward to check that Tutte(:e, y) € #P when x, y are 
integers with x,y > 1. If cc,y are arbitrary integers then the terms in the 
Tutte polynomial vary in sign, and the problem Tutte(:e, y) no longer fits 
the #P framework. In that case, Tutte(x, y) E GapP, where GapP is 
the set of functions / = / + — f~ : T,* — > Z expressible as the difference 
of two #P-functions / + and (Simply partition the terms of the Tutte 
polynomial according to sign, and compute the positive and negative parts 
separately.) 

Finally, since we do not want to restrict ourselves to integer x and y, we 
need to extend the classes #P and GapP a little to encompass computations 
over the rationals. We say that / : X* — > Q is in the class #Pq if f(x) = 
a(x)/b(x), where a, b : X* — ► N, and a € #P and b G FP, where FP is the 
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class of functions computable by polynomial-time algorithms. If x,y > 1, 
then Tutte(x, y) € #Pq, since we may multiply through by suitable powers 
of the denominators of x and y, after which all the terms in the Tutte 
polynomial become integers. 

A randomised approximation scheme is an algorithm for approximately 
computing the value of a function / : X* — ► R. The approximation scheme 
has a parameter e > which specifies the error tolerance. A randomised 
approximation scheme for / is a randomised algorithm that takes as input 
an instance x € X* (e.g., an encoding of a graph G) and an error tolerance 
e > 0, and outputs a number z 6 Q (a random variable of the "coin tosses" 
made by the algorithm) such that, for every instance x, 

Pr [e~ £ f(x) < z < e £ f(x)] > | . (2) 

The randomised approximation scheme is said to be a /uZZy polynomial ran- 
domised approximation scheme, or FPRAS, if it runs in time bounded by a 
polynomial in \x\ and e _1 . Note that the quantity 3/4 in Equation ([2]) could 
be changed to any value in the open interval (|, 1) without changing the set 
of problems that have randomised approximation schemes [13} Lemma 6.1]. 

It is known that every counting problem in #P has a randomised ap- 
proximation scheme whose complexity is not much greater than NP. In 
particular, if / is a counting problem in #P then the bisection technique of 
Valiant and Vazirani [211 Cor 3.6] can be used to construct a randomised 
approximation scheme for / that runs in polynomial time, using an oracle 
for an NP predicate. See |13l Theorem 3.4] or [5l Theorem 1]; also [18] for 
an early result in this direction. 

We will use the notion of approximation-preserving reduction from Dyer, 
Goldberg, Greenhill and Jerrum [5] . Suppose that / and g are functions from 
X* to R. An "approximation-preserving reduction" from / to g gives a way 
to turn an FPRAS for g into an FPRAS for /. An approximation-preserving 
reduction from / to g is a randomised algorithm A for computing / using 
an oracle for g. The algorithm A takes as input a pair (x,e) € X* x (0, 1), 
and satisfies the following three conditions: (i) every oracle call made by A 
is of the form (w, 5), where w £ X* is an instance of g, and < 5 < 
1 is an error bound satisfying <5 -1 < poly(|»|, (ii) the algorithm A 

meets the specification for being a randomised approximation scheme for / 
(as described above) whenever the oracle meets the specification for being 
a randomised approximation scheme for g; and (iii) the run-time of A is 
polynomial in \x\ and s^ 1 . 

If an approximation-preserving reduction from / to g exists we write 
/ <AP 9, an d say that / is AP-reducible to g. Note that if / <ap g and g 
has an FPRAS then / has an FPRAS. (The definition of AP-reduction was 
chosen to make this true). If / <ap g and g <ap / then we say that / and 
g are AP-interreducible, and write / =ap Q- 
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Dyer et al. [5] identified three classes of counting problems that are 
interreducible under approximation-preserving reductions. The first class, 
containing the problems that admit an FPRAS, are trivially AP-interreducible 
since all the work can be embedded into the reduction (which declines to use 
the oracle). The second class (and the last one that we will describe here) 
is the set of problems that are AP-interreducible with #Sat, the problem 
of counting satisfying assignments to a Boolean formula in CNF. Zucker- 
man |24j has shown that #Sat cannot have an FPRAS unless RP = NP. 
The same is obviously true of any problem in ^P to which #Sat is AP- 
reducible. See [5] for details. 

2.3 The Tutte polynomial and #P 

We will study the following computational problem for rationals x and y. 
Name Tutte(x, y). 
Instance A graph G = (V,E). 
Output T(G;x,y). 

Given fixed rationals x and y, Tutte(x, y) is a function from S* to 
Q, mapping an encoding of a graph G to a rational T(G;x,y). It is not 
immediately clear from the definition (pQ) that TuTTE(a;,y) is in #Pq, but 
this is known to be true if x and y are both non-negative. 

In particular, if G is connected, it is known [19] (see also [20 \ Theorem 
1X.65]) that T(G;x,y) can be expressed as T(G;x,y) = Ysr x a ^y b ^ T \ 
where the sum is over spanning trees T of G, and a(T) and b(T) are natural 
numbers which are easily computable from T0 (a>(T) is the number of so- 
called "internally active" edges of T and b(T) is the number of "externally 
active" edges of T — see [20] for details.) 

It is clear from the definition (pQ) that the Tutte polynomial of a graph 
G (which may have several connected components) may be expressed as a 
product of the Tutte polynomials of the components. Thus, for x > and 
y > 0, we have Tutte(x, y) £ #Pq, which means that there is a randomised 
approximation scheme for Tutte(x, y) that runs in polynomial time, using 
an oracle for an NP predicate. 

It is unlikely that Tutte(x, y) is in #Pq for all x and y. In particular, 
Theorem [T7] identifies a region of points (x, y), where y is negative, for which 
even approximating Tutte (x, y) is as hard as #P. 

2.4 Previous work on the complexity of the Tutte polynomial 

Jaeger, Vertigan and Welsh [iO] have completely mapped the complex- 
ity of exactly computing the Tutte polynomial. They have observed that 

1 Indeed, historically, this appears to have been the original definition of the polynomial. 
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Tutte(x, y) is in FP for any point (x,y) on the hyperbola Hi. This can 
be seen from the definition (pQ), since terms involving k(A) cancel. Also, 
Tutte(x, y) is in FP when (x, y) is one of the special points (1, 1), (0, —1), 
(—1,0), and (—1,-1). As noted in Section Q3 T(G; 1, 1) is the number of 
spanning trees of a connected graph G, T(G; 0, —1) is the number of 2-flows 
of G (up to a factor of plus or minus one), and T(G; —1,0) is the number 
of 2-colourings of G (up to an easily computable factor). T(G; —1, —1) has 
an interpretation in terms of the "bicycle space" of G. See \10\ (2.8)]. In- 
triguingly, Jaeger, Vertigan and Welsh managed to show that Tutte(x, y) 
is #P-hard for every other pair of rationals (x,y). They also investigated 
the complexity of evaluating the Tutte polynomial when x and y are real or 
complex numbers, but that is beyond the scope of this paper. 

The only FPRAS for approximating the Tutte polynomial that we know 
of is the Ferromagnetic Ising FPRAS of Jerrum and Sinclair [12] . This gives 
an FPRAS for Tutte(x, y) for every point (x,y) on H2 with y > 1. The 
connection between the Ising model and the Tutte polynomial along the 
hyperbola H2 is elaborated later in the paper — see Equation (f39|) . We 
know of no other FPRASes for approximating the Tutte polynomial for an 
arbitrary input graph G. There is some related work, however, for example, 
Karger [14J gives an FPRAS for non-reliability, which is not the same thing 
as an FPRAS for reliability, but is somewhat related. 

There are also FPRASes known for special cases in which restrictions 
are placed on G. For example, pQ gives an FPRAS for points (x, y) with 
x > 1 and y > 1 for the restricted case in which the input graph G is 
"dense", meaning that the n- vertex graph G has minimum degree Q(n). 
As another example, there is a huge literature on approximately counting 
proper colourings of degree-bounded graphs. 

Several negative results are known for approximating the Tutte polyno- 
mial. First, note that if T(G; x, y) is the number of solutions to an NP- 
complete decision problem, then there can be no FPRAS for Tutte(x, y) 
unless RP = NP. So, for example, if A > 2 is a positive integer, then by 
the chromatic polynomial specialisation mentioned in Section [Jj there is no 
FPRAS for T(G;1 - A,0). 

Jerrum and Sinclair |12} Theorem 14] showed that there is no FPRAS 
for (antiferromagnetic) Ising unless RP = NP. This implies that, unless 
RP = NP, there is no FPRAS for the function whose input is a graph G 
and a point (x, y) on H2 with < y < 1 and whose output is T(G; x, y). 

Welsh [23] extended this result. Specifically, he showed the following, 
assuming RP 7^ NP. 

• Suppose q > 2 is a positive integer. Then there is no FPRAS for 
the function whose input is a graph G and a point (x, y) on H q with 
x < 0, y > and whose output is T(G; x, y). Furthermore, there is no 
FPRAS for the function whose input is a graph G and a point (x,y) 
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on H q with x < 0, y < and whose output is T(G; x, y). 

• There is no FPRAS for the function whose input is a graph G and a 
point (x,y) on #3 with < x < 1 and whose output is T(G;x,y). 

3 Regions of the Tutte plane that do not admit an 
FPRAS unless RP=NP 

The tensor product of matroids was introduced by Brylawski [2]. We define 
it here in the special case of graphs. Let G be a graph, and K another graph 
with a distinguished edge e with endpoints u and u' . The tensor product 
G®K is obtained from G by performing a 2-sum operation with K on each 
edge / of G in turn: Let the endpoints of / be v and v'. Take a copy of K 
and identify vertex u (resp. u 1 ) of K with v (resp. v') of G, and then delete 
edges e and /. (Since G and K are undirected graphs, there are two ways 
of performing the 2-sum. This lack of uniqueness is an artefact of viewing 
a matroid operation in terms of graphs, which have additional structure. 
However, the Tutte polynomial is insensitive to which of the two possible 
identifications is made.) For technical reasons we will assume that e is not 
a bridge of K. In particular, we assume that deleting e does not increase 
the number of connected components of K. 

Let K\e be the graph constructed from K by deleting edge e. Let K/e be 
the graph constructed from K by contracting edge e. Suppose (x,y) 6 Q 2 . 
Let q = {x — l)(y — 1). Define the point (x',y') as follows. 



and 



/ = (l-q)T(K\e;x,y) 

T(K\e;x,y)-(x-l)T(K/e;x,y) {) 

/ _ (1 - q)T(K/e;x,y) 

T(K/e;x,y)-(y-l)T(K\e;x, y y {) 



Then it is known (pUl (4.1)]) that 

T(G; x', y') = L(x, y, K) m M(x, y, K) n ~ K T(G ® K; x, y), (5) 

where n, m and k are (respectively) the number of vertices, edges and con- 
nected components in G and 

L(x,y,K) 



T(K/e;x,y)-(y-l)T(K\e;x, y y 
and 

M( R) = T(K/e;x,y)-(y-l)T(K\e;x,y) 
V ' y ' ; T(K\e;x,y)-(x-l)T(K/e;x,y)- 
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Suppose that the denominators of ([3j) and are non-zero. In this case, 
the point (a/, y') is well-defined and we say that (x, y) is shifted to the point 
(x',y') by K. In this case, L(x,y,K) and M(x,y,K) are also well-defined, 
so Equation (JSJ) gives us the reduction Tutte(x', y') <ap Tutte(z, y). 

We will not prove Equation ([5]) since the equation can be found elsewhere 
(eg., pUl (4.1)]) but, for completeness, we will derive similar identities that 
we will use in Section 01 We are particularly interested in two special cases. 
The case in which K is a cycle on k + 1 vertices is known as a k-stretch in the 
literature and the case in which K is a two- vertex graph with k + 1 parallel 
edges is known as a k-thickening. Informally, a /c-stretch of G replaces each 
edge of G by a path of length k, while a A;-thickening replaces each edge by 
a bundle of k parallel edges. Specifically, 



Observe that q = (x — l)(y — 1) is an invariant for stretches and thick- 
enings, and indeed for shifts in general. It is this limitation that gives the 
hyperbolas H q a special place in the complexity theory of the Tutte poly- 
nomial. All shifts preserve q = (x — l)(y — 1) but not all AP-reductions do. 
In particular, the construction in [TU1 (5.12)] (taking p = 1), based on an 
idea of Linial pj)], gives the reduction Tutte(x, 0) <ap Tutte(x — 1, 0) for 



We shall make frequent use of the fact that shifts may be composed. 
Lemma 1. The relation "shifts to" is transitive. 

Proof. Suppose K\ is a graph that implements the shift (x,y) — > (x',y') 
and K 2 is the graph (with distinguished edge e) that implements (x',y') — > 
(x",y"). Let K be the graph obtained from K2 by performing a 2-sum 
with K\ along every edge of K2 except e; let e remain the distinguished 
edge of K. We claim that K implements the shift (x,y) — > (x",y"). Since 
G®K = (G <S> K2) <S> K\, for any G, this ought to be true, but we can verify 
the claim by direct calculation. 

Evaluating the rhs of ([3]), with K = K: 



(l-g)T(K\e;x,y) 
T(K\e;x,y)-(x-l)T(K/e;x,y) 

= (l-q)T{{K2\e)®K l] x : y) 

T((K 2 \ e) ® Ki);x, y) - (x - l)T((K 2 /e) ®K 1]X , y) 

^ (l-q)T(K2\e;x>,y>) 

T(K 2 \ e; x>, y>) - (x - l)M(x, y, K x )T{K 2 /e; x>, y>) 

(l-q)T(K2\e;x',y') 
T(K 2 \ e; x', y>) - (x> - l)T(K 2 /e; x', y>) 




(x k ,q/(x k — 1) + l) for a /c-stretch; 
{?l/{y k ~ 1) + l>y fc ) f° r a ^-thickening. 



(6) 



x/1. 



= X' 



,11 



(9) 
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Here, © uses ©, and © the fact that (x' - 1) = (x - l)M(x,y,Ki). 
Equality ([7J) follows from ([5]), noting that K^/e has the same number of 
edges as K% \e, but one fewer vertex. A similar calculation holds for y" . □ 

Shifts play a key role in the classical study of the complexity of exact 
computation of the Tutte polynomial |10j . and the same is true in the current 
investigation. Our keys tools are the following. 

Theorem 2. Suppose (x,y) G Q 2 satisfies q = (x — l)(y — 1) ^ {0,1,2}. 
Suppose also that it is possible to shift the point (x,y) to the point (x',y') 
with y' [-1, 1], and to (x", y") with y" G (-1, 1). Then there is no FPRAS 
for Tutte (x, y) unless RP = NP. 

Theorem 3. Suppose (x,y) G Q 2 satisfies q = (x — l)(y — 1) ^ {0,1,2}. 
Suppose also that it is possible to shift the point (x, y) to the point (x', y') 
with x' $l [—1, 1], and to (x" , y") with x" G (—1, 1). Then there is no FPRAS 
for Tutte(x, y) unless RP = NP. 

Since the notion of "shift" is defined for any class of matroids closed 
under tensor product, it should be possible to frame statements similar to 
Theorems [2] and [3] for classes of matroids other than graphic. Although the 
two theorems are dual to one another in the matroid theoretic sense, they 
are not equivalent, since the class of graphic matroids is not closed under 
duality. 

The proofs of Theorem [2] and [3] appear in Section 01 
3.1 Two halfplanes 

Corollary 4. Suppose (x,y) is a point lying in the open half-plane x < — 1 
but not on the hyperbolas Hq or H\. Under the assumption RP ^ NP there 
is no FPRAS for Tutte(x, y) . 

Proof. Let (x,y) £ I 2 be a point not on Hq or H\ that satisfies x < — 1. 
At the outset, we'll assume further that (x,y) £ Hi and that y ^ — 1. 
There are three cases, depending on y. First assume y > 1, and observe 
that q = {x — l)(y — 1) < 0. Using a /c-stretch, we may shift the point 
(x,y) to the the point (x",y") = (x k ,q/(x k - 1) + l). Now y" G (-1,1) 
for all sufficiently large even k so Theorem [2] applies. (The trivial shift, 
taking (x,y) to itself, provides the point (x',y') with y' ^ [—1, 1].) A similar 
argument, but setting k to be large and odd deals with the situation y < —1. 
Finally, when y G ( — 1,1), a 2-stretch shifts (x,y) to the point (x',y') = 
(x 2 , q/(x 2 - 1) + 1) = (x 2 , (y - l)/(x + 1) + l) , with y' > 1. 

The additional condition y ^ — 1 may be removed by noting that a 3- 
stretch shifts (x,— 1) to a point (x" \y") = (x 3 , 1 — 2/(x 2 + x + 1)) with 
x" < —1 and y" G (— 1,+1), and we have already seen how to deal with 
such a point. 
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Finally, suppose q = 2. Like Welsh [23] we will show hardness using an 
argument of Jerrum and Sinclair |12l (Theorem 14)]. Suppose that G has n 
vertices and m edges and that x' and y' satisfy {x' — l)(y' — 1) = 2. Jerrum 
and Sinclair note that 

m 

T(G; x\ y') = (y' - l) n {x> - 1)"< £ ) ]T N r {yT~ r 

r=0 

where N r is the number of functions a : V — > {—1,1} with r bichromatic 
edges. The reader can verify this claim by looking ahead to Equations (I10p 
and (|39p . Thus, if G has a cut of size b then 

T(G;x',y') > (y' - l) n {x> - l)~^ E \yT- b - 

Otherwise, 

T{G;x',y') < (y 1 - l) n (x' - i)-^m 2 n ( y ') m - b +\ 

Now consider a point (x, y) on H2 with x < —1. Note that y 6 (0, 1). Let 
A; be a positive integer with y k < 2~ 2n and let y' = y k . Let x' = 2/(y k — 1) + 1. 
If we had an FPRAS for Tutte(x, y), we could estimate T(G;x',y') by k- 
thickening. Thus, we could determine whether or not G has a cut of size b, 
giving RP = NP. 

□ 

Corollary 5. Suppose (x,y) is a point lying in the open half-plane y < —1 
but not on the hyperbolas Hq, Hi or H2. Under the assumption RP 7^ NP 
there is no FPRAS for TuTTE(j;,y). 

Proof. Dual to the proof of Corollary 2] (but without the extra argument for 
q = 2). □ 

Corollaries H] and exclude the hyperbola q = 0. Nevertheless, the 
arguments of Theorem [2] can be extended to handle the portion of this 
(degenerate) hyperbola in which y < — 1. Specifically, in Section U] we prove 
the following. 

Lemma 6. Suppose (x, y) is a point with x = 1 and y < — 1. Then there is 
no FPRAS for Tutte(x, y) unless RP = NP. 

We do not know whether the arguments of Theorem [3] can be similarly 
extended to q = 0. 

The hyperbola H2 is excluded from Theorems [2] and [3] and a separate 
argument (following Welsh) was used to include H2 within the scope of 
Corollary H] which applies to the region x < —1. We do not know of a 
similar argument which applies to H2 in the region y < — 1 and indeed this 
hyperbola seems to have a special status in the region y < —1, as Lemma [7] 
shows. Consider the following computational problem. 
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Name ^Perfect Matchings. 
Instance A graph G. 

Output The number of perfect matchings in G. 

^Perfect Matchings is #P-complete, but it is not know whether it has 
an FPRAS. In Section H] we prove the following. 

Lemma 7. Suppose (x, y) is a point on the hyperbola Hi with y < — 1. Then 
TuTTE(x,y) =ap #Perfect Matchings. 

Remark: For convenience, we allow the graph G in the definition of 
^Perfect Matchings to have loops and/or multiple edges. This is with- 
out loss of generality, since the perfect matchings of a graph G are in one-to- 
one correspondence with the perfect matchings of the 3-stretch of G (which 
has no loops or multiple edges). 



3.2 The Vicinity of the Origin 

In this section, we consider the region given by \x\ < 1 and \y\ < 1. We 
have already seen (in the proof of Corollary 0]) that, unless RP = NP, 
there is no FPRAS for Tutte(;c, y) for any point (x, y) on the hyperbola 
(x — l)(y — 1) = 2 in this region. The following lemmas give additional 
regions that do not admit an FPRAS unless RP = NP. 

Lemma 8. Suppose (x,y) is a point with \x\ < 1, \y\ < 1 and y < —1 — 
1x that does not lie on the hyperbola Hi. Then there is no FPRAS for 
Tutte(x, y) unless RP = NP. 

Proof. Using Equation ©, a 2-stretch shifts (x,y) to (x',y') with 

, _ (a - !)(?/ - 1) _ y + x 
x z — 1 X + 1 

Now if q = {x — l)(y — 1) ^ {0, 1, 2}, the lemma follows from Theorem^ As 
we noted above, the result is already known for q = 2. Also, Ho is outside 
the scope of the lemma. □ 

Lemma 9. Suppose (x,y) is a point with \x\ < 1, \y\ < 1 and x < — 1 — 2y 
that does not satisfy {x — l)(y — 1) = 1. Then there is no FPRAS for 
Tutte(x, y) unless RP = NP. 

Proof. As in the proof of Lemma [8] we can use a 2-thickening together with 
Theorem[3]to obtain the result when except for q = 0, q = 1 and q = 2. The 
result is known for q = 2 and the cases q = and q = 1 are excluded from 
the lemma. □ 
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Lemmas [8] and [9] give two intractable open triangles in the vicinity of the 
origin. The following lemmas extend intractability to the boundaries. The 
value 0.29 in the statement of the lemmas has no special significance. We 
do not know whether the entire boundary is intractable, but the value 0.29 
is not best possible — it was chosen because it yields a simple proof. 

Lemma 10. Suppose (x,y) is a point with x = — 1 and —\<y< 0.29, 
excluding the special point (x,y) = (—1,0). Then there is no FPRAS for 
Tutte(x, y) unless RP = NP. 

Proof. A 2-thickening of (— 1, y) gives the point (x', y 1 ) = ((y — l)/(y+l), y 2 ). 
If y < then x' < -1 so the result follows from Corollary 0] since (x',y') is 
not on Hq or Hi. Now if < y < 1 then x' G (—1, 0) so \x'\ < 1 and \y'\ < 1. 
Now note that if < y < 0.29 then y' < —1 — 2x' so the result follows from 
Lemma El □ 

Lemma 11. Suppose (x,y) is a point with y = — 1 and — 1 < x < 0.29, 
excluding the special point (x,y) = (0,-1). Then there is no FPRAS for 
Tutte(x, y) unless RP = NP. 

Proof. Dual to the proof of Lemma [TU1 □ 

The intractable triangles from Lemma [8] and [9] certainly do not cover all 
intractable points in the vicinity of the origin. Possibly the whole of the 
region \x\, \y\ < 1 is intractable (apart from Hi and the special points). 

Here is a lemma which adds a little bit to our knowledge in the region. 
For example, it includes the point (x,y) = (—0.23, —0.23) which has q > 1.5 
but is not covered by Lemma [8] or [9j 

Lemma 12. Suppose (x,y) is a point with \x\ < 1 and \y\ < 1 and (x — 
l)(y — 1) = q > 1.5 (excluding the special points (—1,-1), (—1,0) and 
(0, -1) ). Then there is no FPRAS for Tutte(x, y) unless RP = NP. 

Proof. First, note that neither x nor y is 1, since that would make q = 0. 
Also, we don't have (x,y) = (—1, —1) since that is a special point. Suppose 
x = —1. Then y > —1. The restriction on q implies y < 0.25, so the result 
follows from Lemma [TU1 Similarly, the case y = — 1 follows from Lemma [TT1 
So suppose \x\ < 1 and \y\ < 1. 

If q > 2 then the result follows from Theorem O Do a 2-thickening 
(Equation ([6])) to shift to the point 



Note that y 2 — 1 G ( — 1,0) since \y\ < 1 so q/(y 2 — 1) < — q < —2 since 
q > 2. So x' < —1. Then apply Theorem [3j The case q = 2 is known, as 
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noted at the beginning of the section. Suppose 3/2 < q < 2. For a large 
even integer k, do a /c-thickening to shift (x, y) to the point 



Choose k so that < y k < (2 — q)/2 (this is possible since q < 2). Consider 
x\ — 1 = q/(y k — 1). Note that this is in the interval (—2, — q). Now do a 
2-stretch to shift (xi,yi) to the point (x',y r ) = (x\,q/{x\ — 1) + 1). Note 
that 



where the upper bound of —2 follows from the bounds that we derived on 
x\ and q > 3/2. Now use Theorem [2j □ 

The lemma could certainly be improved. For example, consider the point 
(x,y) = (—0.2,0) with q = 1.2. An alternating sequence of 14 2-stretches 
and 2-thickenings shifts this point to a point (x',y') ~ (—103.1,0.99) so 
(x,y) had no FPRAS (unless RP = NP) by Theorem El 

4 The reductions 

4.1 The Multivariate Formulation of the Tutte Polynomial 

It is convenient for us to use the multivariate formulation of the Tutte poly- 
nomial, also known as the random cluster model [22} I17j . For a graph 
G = (V, E) with edge weights w : E — > Q and q € Q, define the multi- 
variate Tutte polynomial of G to be Z(G;q,w) = ^2ace w (^) < 1 k( ' A ^ where 
w(A) = JleeA w ( e )i an d is the number of connected components in the 
graph (V,A). 

Suppose (x,y) G Q 2 and q = (x — l)(y — 1). For a graph G = (V,E), 
let w : E — > Q be the constant function which maps every edge to the value 
y — 1. Then (see, for example pH (2.26)]) 

T(G; x, y) = (y- iy n (x - 1)-^Z(G; q, w). (10) 

So Z is a generalisation of T that allows different weights to be assigned 
independently to different edges. For rationals q and 7, let MultiTutte(q, 7) 
be the following problem. 

Name MultiTutte((/, 7) . 

Instance A graph G = (V, E) with edge labelling w where w is the constant 
function mapping every edge to the value 7. 

Output Z(G;q,w). 
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Suppose (g, 7) € Q 2 . Equation (fTUj) gives us the reduction 



MULTlTUTTE(g, 7) <AP TUTTE + 1,7+ 1 

If (», y) G Q 2 and neither x nor y is 1 then Equation (fTOj) gives us the 
corresponding reduction 

Tutte(x, y) <ap MultiTutte((x - l)(y - l),y — 1). 

Not surprisingly, the notion of a shift from S|3] may be re-expressed in 
terms of the new parameters. Doing so has the advantage of allowing us to 
apply shifts to individual edges of a graph, as opposed to the whole graph. 
This idea is explored in [17\ §4.6]. We derive the equations that we need 
here in order to fix the notation and explore the concepts. As in JJJ let K 
be a graph with distinguished edge e, and suppose that K has constant edge 
weight a S Q. Define 



/ _ qZ(K/e;q,a) - Z(K\e;q,a) 
a Z(K\e;q,a)-Z(K/e;q,a) [ ] 



and 

N(q,a,K) 



q(q-l) 



Z(K \ e; q, a) - Z(K/e; q, a) 



Let G = (V, E) be a graph, / € E and w' : E — > Q a weighting such that 
w'(f) = a' . Denote by Gf the 2-sum of G and K along /. Let w be the 
weight function on Gf that gives every edge of K weight a and inherits the 
remaining weights from w' . We will show below that 

Z(G; q, w') = N(q, a, K) Z(G f ;q, w). (12) 

One way to capture (|12p informally is to say that a single edge of weight a' 
may be simulated by a subgraph K whose edges have weight a. 

Suppose that the denominator of (fTT|) is non-zero. In this case, the point 
(q, a') is well-defined and we say that (q, a) is shifted to the point (q, a') 
by K. In this case, N(q, a, K) is also well-defined, so Equation (fT2"j) gives us 
an efficient algorithm for approximating Z(G; q, w') by using an subroutine 
for computing Z(Gf\q,w). 

For the derivation of (fTTI) and (fT2|) . let v and v' be the endpoints of /. 
Let S be the set of subsets of E — {/} which connect v and v' and let T be 
the set of all other subsets of E — {/}. Then Z(G; q, w') = Z$ + Zt, where 

Z S = ^2w(A')q k ( A '\l+a'), 

A'es 

and 

Z T = £ W (A')q k ^ (l + i 
A'eT 
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Similarly, Z(Gf,q,w) = Zf t s + ^/,T> where 



Z ft a = ™(A')q k(A,) 



A'es 

and 



z f , T =Y / ^ A ')^ A ' ) Z{KX 2 q,a) 



A'eT 

Now the equation for a' comes from the following argument. Suppose we 
could define a' such that 

^ = f^, (13) 

and N(q, a, K) = Zs/Zf g. Then (|12p would hold as desired. Now note that 
(fTTl) entails (flTfl) . 

The shifts that we have defined here are consistent with the usage in £j3l 
In particular, suppose that (x, y) is shifted to the point (x' , y') by a graph K 
with distinguished edge e. As long as e is not a bridge of K then taking 
a = y — 1 and a! = y' — \ and g = (x — l)(y — 1) we find (from Equations (jj]) 
and ([10]) and (fTTI) ) that the same graph K shifts (g, a) to (q,a'). 

Thus, the equation describing stretching and thickening, Equation ©, 
can be translated as follows. (See, for example, pU (4.20), (4.26)]) 

q \ k 



— = hi — 1, for a /c-stretch; , , 

a' \a J (14) 

a' = (a + l) k — 1, for a A;-thickening. 

We now generalise the computational problem MultiTutte((7, 7) de- 
fined earlier. For rationals q, a\, . . . , a^, MultiTutte((7; ati,..., a^) is the 
problem: 

Name MultiTutte((/; a\, . . . , «&). 

Instance A graph G = (V, E) with edge labelling w : E — > {qi, . . . , a^}. 
Output Z(G;q,w). 



4.2 Proof of Theorem G3 

The decision problem Minimum 3- way Cut is: 

Instance A simple graph G = (V, E) with three distinguished vertices 
("terminals") t\,t2,t^ € V, and an integer bound b. 

Output Is there a set of at most b edges whose removal from G disconnects 
ti from tj for every i, j £ {1, 2, 3} with i 7^ jl 

It was shown to be NP-complete by Dahlhaus et al. [3]. 
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Lemma 13. Suppose q G Q \ {0,1,2}, and that a\,a2 € Q satisfy a\ ^ 
[—2,0] anda2 G (—2,0). T/ien i/iere is no FPRAS for MultiTutte(<;; ai, 02) 
un/ess RP = NP. 

Proof. Suppose G = (V, E,ti,t2,ts) is an arbitrary instance of Min 3-WAY 
Cut. Without loss of generality assume G is connected, and for conve- 
nience let n = \V\ and m = \E\. Our ultimate goal is to construct an 
instance (G',w') of MultiTutte(o/; 0:1,02) such that Z(G';q,w') is a close 
approximation to the number of minimum 3-way cuts in G. The size of a 
minimum cut will be a by-product of the of the reduction. 

As an intermediate goal, we'll construct a weighted graph (G = (V, E),w) 
such that Z{G; q, w) is a close approximation to the number of minimum 3- 
way cuts in G where w : E — > {0\,02} for some conveniently-chosen values 
0i and 02- The final step of the proof will be to relate these convenient 
values to the specified ones, namely a\ and 02- We will require 0\ to be 
sufficiently large, in particular, let q = max(|oJ, 1). Let M = 8 x 2 m q n . We 
will require 

0i > M. (15) 

We will also require 02 to be sufficiently close to —1. In particular, we will 
choose a small value 5 (see Equations (f2Tj) and (f22|) . depending on m, q and 
n. We will require [1 + 02] < S. 

The construction of (G = (V, E),w) is very direct: V = V, E = EUT, 
where T = {{h,t 2 }, {^2,^3}, ^3}}, and 



w{e) 



01, ifeeE; 

02, otherwise. 



Now, letting 

Ai\2,3 = {A C E : h ^Ah and h </> A h and t 2 ~a t 3 }, 

etc, where ~a denotes the binary relation "is connected to" in the graph 
(V, A), we may express the multivariate Tutte polynomial of G as 

Z(G; q, w) = Zi| 2 | 3 + A|2,3 + S 2\l,3 + S 3\l,2 + £1,2,3, (16) 

where, e.g., 

The overview of the proof is as follows: we show that for 0\ sufficiently 
large and 02 sufficiently close to —1, the last four terms on the r.h.s. of (|16|) 
are negligible in comparison with the first, and that the first term, -^"1 1 2 1 3 , 
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counts minimum 3-way cuts in G (approximately, and up to an easily com- 
putable factor). Up to symmetry there are three essentially distinct terms 
in (|16p . and we consider them in turn. First, 



1^1 



2,3 



52 ^w(AU B)q^ AuB ^ 
AeAi,2,3 bct 



52 w(A)q« A \l + p 2 ) 
<(2p 1 ) m \q\q n - 1 d 3 , 



Here we have used 1 < n{A) < n. Next, 



1 2,31 



52 52w(AuB)q^ AuB ^ 
AeA ll2 , 3 BCT 



52 w{A)q^- l {q + 2(5 2 +Pl){l+f5 2 ) 
AeAi\ 2 , 3 

< {2^) m \q\q n ~ 2 (g + 4 + 4)5 

< 9(2/3 1 ) m |g|g n - 1 ,5. 



Here we used \f3 2 \ < 2 and 2 < k(A) < n. Last, 



where 
Note that 



i| 2 |3= E E^ U ^ K(AUB) 

Ae-4i|2| 3 BCT 

= £ wi A)q^- 2 ( q 2 +3p 2 q + 3f3 2 2 +Pl) 
AeA 1}2 \3 

AeA ll2 \3 



C{(3 2 ) = (q - l)(g - 2) + 3(g - 1)(1 + (3 2 ) + (1 + (3 2 f 



\C(f3 2 )-(q-l)(q-2)\<3\q-l\5 + 5 3 , 



(17) 



(18) 



(19) 



The crucial fact is that C{j3 2 ) remains bounded away from as S — > 
(and hence (3 2 — > —1), provided (as we are assuming) q ^ {1,2}, whereas 
expressions (fTT|) and (JTHJ) tend to as 5 — > 0. 

(i) 

Now denote by -4^|2|3 the set of all subsets in *4i|2|3 of size i. Let c be 
the size of a minimum 3-way cut in G, and N be the number of such cuts. 
Then 
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27- 



1|2|3 



m—c—1 



AeA 



Am-c) i=0 
**1|2|3 ^ lt -^l|2|3 

m—c—1 



i=0 AeA (i) 



Then 



is equal to 



A|2|2 



N 



' m—c—1 



' =0 ^2,3 



k(A)-2 



(/3T C ?) 



so crudely upper-bounding the absolute value of the right-hand-side, we get 



1 1 2 1 3 



c(/3 2 )/?r- c <? 



iV 



< 



2 m q n 



Now set <5 to satisfy both 



and 



icroi > ||C(-i)] 

5< 



(20) 



(21) 



(22) 



448(2/3i)™g" 

Now (JED ensures that the r.h.s. of flU) is at most ± while (fTT]) . (fTHI). (fT9 



21]) and ([22]) ensure 

A|2|3 - Z(G;q,w) 



A|2,3 + ^2|1,3 + ^3\1,2 + A, 2,3 



< 



Combining this inequality with (|20|) . the bottom line is 

Z(G;q,w) 



N 



1 

< -. 

- 4 



(23) 



If we knew Z(G;q,w), we could determine c — it is the unique integer 
such that (|23|) provides an estimate for N that lies in the range [1, 2 m ]. The 
value of c is unique since j3\ > M > 2 m .) 

In fact, we do not need an exact value of Z(G;q,w) — an approxi- 
mate value will do. In particular, an FPRAS for Z(G; q, w) would give a 
randomised polynomial-time algorithm for computing c, which would show 
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RP = NP. For details about approximation accuracy, see [I], especially the 
final three paragraphs of the proof of Theorem 3. 

Finally we need to relate our weights fli, [3 2 to the given ones ai, ai- Let 
positive integers ki,k 2 satisfy [a\ + l) kl — 1 > M and Ua 2 + l) fc2 | < $■ 
Let K\ be a 2- vertex graph with k\ + 1 parallel edges, each of weight a±. 
Recall that taking a 2-sum with K\ implements a k\ thickening. Let K 2 be 
a 2- vertex graph with k 2 + 1 parallel edges, each of weight a 2 . Let G' be 
the graph derived from G by taking the 2-sum of each weight 0\ edge with 
K\ and taking the 2-sum of each weight (5 2 edge with K 2 . Call the resulting 
graph G' and its weighting w' . By repeated application of (fT2l) . Z'(G ! ; g, = 
N(q,a 1 ,K 1 ) m N(q,a 2 ,K 2 ) 3 Z{G';q,w'). By setting A = fa + l) fc * - 1, for 
i = 1, 2, we satisfy /?i > M and |/?2 + 1| < S, as required by our reduction. 
(This is by (|14p and the definitions of k± and A^.) Finally observe that 
k\ = 0(m) and A;2 = 0(m 2 ), so the size of G' is polynomially bounded. 

Thus an FPRAS for MuLTlTuTTE(g; a\, a 2 ) would yield a polynomial- 
time randomised algorithm for computing the size of a minimum 3- way cut, 
which would entail RP = NP. □ 

Using Lemma [T3"l we can now prove Theorem [2l 
Theorem H Suppose (x,y) E Q 2 satisfies q = (x - l)(y - 1) ^ {0,1,2}. 
Suppose also that it is possible to shift the point (x, y) to the point (x 1 , y') 
with y' (£ [-1, 1], and to (x", y") with y" G (-1, 1). Then there is no FPRAS 
for Tutte(x, y) unless RP = NP. 

Proof. Let a = y — 1 and a\ = y' — 1 and a 2 = y" — 1. Note that a\ [—2, 0] 
and a 2 £ (—2,0). Let (K',e') be a graph that shifts (x,y) to (x',y r ) and 
note that (K',e') also shifts (g, a) to (q, ot\). Similarly, suppose (K",e") 
shifts (x,y) to (x",y") and therefore shifts (g, a) to (q,a 2 ). 

Suppose (G,w) is an instance of MuLTiTuTTE(g; ai, a 2 ) with mi edges 
with weight a\ and m2 edges with weight a 2 . Denote by G the graph derived 
from G by taking a 2-sum with (K', e') along every edge with weight «i and 
taking a 2-sum with (iT", e") along every edge with weight a 2 . Let w be the 
constant weight function which assigns weight a to every edge in G. 

Then by repeated use of Equation (fT2|) . 

Z(G; q, w) = N(q, a, K x ) m ^N(q, a, K 2 ) m2 Z(G; q, w). 
Thus by Equation ([ID]) . 

Z(G; g, «;) = iV(g, a, Ki) mi iV(g, a, K 2 ) m2 (y - l) n (* - 1) K T(G; x, y), 

where n is the number of vertices in G, and k is the number of connected 
components in G. 

Thus an FPRAS for Tutte(x, y) would yield an FPRAS for the problem 
MULTlTuTTE(g; ai,a 2 ), contrary to Lemma [T3"l □ 
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4.3 Extending to q = 

Formally, the multivariate Tutte polynomial Z{G\ q, w) = Y^ace w(A)q K ^ 
is not very interesting at q = because k(A) > 1 so Z(G;q,w) = 0. 
Sokal [T7] treats the q = case as a limit, but for the purpose of ap- 
proximation complexity it is more convenient to work with the polynomial 
Z(G;q,w)q~ K ( E \ We will focus on the case in which G is connected, so we 
define R(G;q,w) = Z(G;q,w)q~ 1 . Note that 

R(G;0,w)= Yl W ( A )' ( 24 ) 

ACE:k(A)=1 

This is the reliability polynomial, and corresponds to the x = 1 component 
of the hyperbola Hq . 

We can express shifts in terms of R{G; q, w). For example, Equation (JT2J) 
does not tell us anything useful about stretching for q = (due to cancel- 
lation) but the same reasoning that we used to derive (jllj) and (]12p gives 
us the following version of these equations for the case in which Gf is a 
A;-stretch (specifically, Gf is the 2-sum of G and a cycle on k + 1 vertices 
along edge /): 

a' = l (25) 

and 

R(G; 0, w') = 1 -L- I R(G f ; 0, w). (26) 

As in the general case, we assume w' is a weight function on G with u>(/) = 
a' and that w inherits its weights from w' except that the new edges in 
the stretch are given weight a. The derivation of (|25|) and (|26p follows the 
derivation of (fTT|) and (fl2l) . Specifically, let S (respectively, T) be the set 
of all subsets A' (1 E — {/} with k(A') = 1 (respectively, k(A') = 2 and 
k{A' U {/}) = 1). Then R(G; 0, u/) = R s + R T where 

^5 = £ w(A')(l + a'), (27) 
A'CS 

i? r = w(A')a'; (28) 

A'CT 

and R(Gf; 0, it/) = -R/,5 + Rf,T> where 

R f S = Y w(A')(a k + ka*- 1 ), (29) 

A'CS 

Rf t T = Yl w ( A ') ak - ( 3 °) 
A'CT 

Similarly, for the case in which Gf is a /c-thickening, we get 

R(G;0,w') = R(G f ;0,w), (31) 
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with a! as in Equation fll4f) 

Now let ZeroMultiTutte(«i, . . . , a fc ) be the following problem. 

Name ZeroMultiTutte^i, . . . 

Instance A connected graph G = (V, E) with edge labelling w : E —* 
{a\, . . . ,a k }. 

Output R(G;0,w). 

An examination of the proof of Lemma [13] gives the following lemma. 

Lemma 14. Suppose that Q.\,a.2 £ Q satisfy ot\ ^ [—2,0] and ct2 £ (—2,0). 
Then there is no FPRAS for ZeroMultiTutte(q!i, a?) unless RP = NP. 

The proof of Lemma [14] follows that of Lemma [T3l By analogy to Equa- 
tion (fTB]) we may express R(G;0,w) as a sum of terms of the form 1 2 1 3 - 
Then 

£i|2|3= w(A)([3l + 3(3 2 2 ), 

Ae^ 1 | 2 | 3 :«(A)=3 

and the other terms all have factors of 5. By analogy to Equation (|20p we 
get 

£ll2l3 2 m 

.V • — . (32) 
Pi 



Using Lemma dU we can now prove Lemma [6] 
Lemma [6], Suppose (x,y) is a point with x = 1 and y < — 1. Then there 
is no FPRAS for Tutte(x, y) unless RP = NP. 

Proof. Let be a point with x = 1 and y < —1. Let a = y — 1 and 

q = 0. Note that a ^ [—2,0]. Let = \_—a\ and let a2 = a/k. Note 
that «2 € (—2,0), and by Equation (f26l) . a /c-stretch shifts (g = 0, a) to 
(<? = 0,02)- Suppose (G,w) is an instance of ZeroMultiTutte(«, a 2 ) 
with 7712 edges with weight «2- Denote by G the graph derived from G by 
applying a fc-stretch to each of these 7712 edges. Let w be the constant weight 
function which assigns weight a to every edge in G. Then by repeated use 
of Equation (|2"o) . 



ma-A).,v) [jJ^t) ma-Mw). 



1112 



Using Equation 



V 7 ACE:k(^) = 1 



23 



where E is the edge set of G, which is connected since G is. Thus, by the 
definition of the Tutte polynomial (pQ), 



where re is the number of vertices of G. So an FPRAS for Tutte(x, y) 
would enable us to approximate R(G;0,w), contrary to Lemma [TU 



4.4 Proof of Theorem [3] 

The following is dual to Lemma [T3l 

Lemma 15. Suppose q E Q \ {0,1, 2}, and that a.\,a,2 £ Q — {0} sat- 
isfy Q/ a i £ [~ 2,0] and q/a2 € (—2,0). Then there is no FPRAS for 
MultiTutte(</; «i, 02) unless RP = NP. 

Proof. We reuse the construction that Frederickson and Ja' Ja' designed in 
order to prove that Connected Bridge-connectivity Augmentation 
(CBRA) is NP-complete (HI Thm 2], though we'll change the edge weights 
to suit our purpose. For convenience the construction will be repeated here. 
We start with an instance of the 3-d Matching Problem: W, X and Y are 
disjoint re-element sets, and MCWxXxYa set of triples. We want 
to know how many "3-d matchings" there are in M. A 3-d matching is a 
subset M' CMofn triples such that every element of W U X U Y is included 
in some triple in M' . For convenience, we'll enumerate the elements of the 
ground set W = {wi, . . . ,w n }, X = {xi, . . .,x n }, and Y = {yi, . . . ,y„}H 

Our ultimate goal is to construct an instance (G' , w') of MultiTutte((/; at\, 
such that Z(G'\ q, w') is determined, to a high degree of accuracy, by the 
number of solutions to the instance of #3-D Matching. In particular, us- 
ing an estimate of Z(G'; q, w'), we'll be able to decide, with high probability, 
whether the number of solutions to the matching instance is zero or strictly 
positive. As an intermediate goal, just as in the proof of Lemma [13| we'll 
construct a weighted graph {G = (V,E),w) that has the desired properties, 
as described above, except that w : V — » {/3i,/?2}> where Pi and fa are set 
to convenient non-zero values. The final step of the proof will be to relate 
these convenient values to the specified ones, namely a\ and a%. The re- 
quirements on 0i and are similar to the ones that we used in the proof of 
Lemma[13l In particular, we will require, for a small e < 1, that \(}\/q\ < e 
(so the absolute value of q/Pi is big). We will also require for a small 5 < 5 
that |1 + qjPi\ < 5 (so P2 is close to —q). We will require e and 5 to be 
sufficiently small — the exact requirements will be given later. 

2 We'll stick, as far as possible, to the notation of [8], though occasional changes are 
needed to avoid clashes. 




□ 
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w 



X Y 




Figure 2: The construction of the graph G in the proof of Lemma [151 The 
edges relating to just one generic triple (wi,Xj,y k ) G M are shown. "Link 
edges" are dashed. 



The vertex set of G (refer to Figure [2]) is 

V = {r}UWUXUYU {a ijk , a ljk : (wi, Xj ,y k ) G M}, 

and the edge set E = T U L where 

T = {{r,Wi},{r,Xi},{r,yi} : 1 < i < n] 

U {{wi,a ijk },{wi,a ijk } : (wi,Xj,y k ) £ M) 
is the set of "tree edges" and 

L= {{xj,a ijk },{a ijk ,a ijk },{a ijk ,y k } : (wi,Xj,y k ) G M). 

the "link edges" . Observe that (V, T) is a tree, and that edges in L join 
leaves in the tree. For e£E, assign weight w(e) = /?2 if e G T is a tree edge, 
and w(e) = Pi if e G L is a link edge. 

We're interested in evaluating Z(G;q,w): 



Z(G;q,w)= J2w(A)q K W 
ACE 

BCLCCT 

= ^/ i (S,/3 2 )/?l B| , (33) 



where 

d\C\ k(BUC) 



h(B;f3 2 ) = E 4 '^ {BUC) = Z(G\B/B;q,w) (34) 
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Here, G \ B/B denotes the graph obtained from G by deleting edges in 
B = L\B and contracting edges in B. 

Let m = \M\, and note that \V\ = 3n + 2m + 1, \T\ = 3n + 2m and 
\L\ = 3m. Our calculation of Z(G;q,w) is greatly simplified if we take (5 2 
to be exactly —q, rather than merely a close approximation. So let's first 
determine, as a function of 5, the absolute error we would introduce by 
replacing (3 2 by —q. Denote by w : E — > Q the weight function 



w(e) 




if e £ T; 
otherwise. 



We wish to estimate the absolute error \Z(G;q,w) — Z(G; q,w)\. Set q = 
max{\q\, 1}; then either |g|/|/?2| > 1, in which case \(3 2 \ < g or |g|/|/^| < 1. 
In this case, since |1 + q//3 2 \ < 1/2, |1 + q/(3 2 \ = 1 - \q\/\l3 2 \ < 1/2, so 
| < 2|g|. We conclude that, in either case, \q\, \f3 2 \ < 2g. Furthermore, for 
all z > 1, we have 

t-i 

01 - (-QT = (ft + q^fti-q)'- 1 -* < mf- 1 ^ = i(2q)% 
3=0 

since \{h + q\ < Iftl |1 + q/fh\ < 2q5. Expanding h(B,/3 2 ) and h(B,-q) 
according to (|34p. and comparing term- by-term, we find that 

\h(B, 02) - h(B, -q) | < 2l T l \T\ (2q) ^Sq l y l 
< |T|(2g) |y|+|T| 5 
= (3n + 2m)(2g) 6n+4m+1 (5. 

So from (f3"3"|) . recalling |/3i| < |g[e < g, 



|Z(G;g,u;)-Z(G;g,™)| < 2l L lgl L l(3n + 2m)(2qf n+4m+1 5 

< (3n + 2m)(2g) 6n+7m+1 <5. (35) 

We'll chose 5 later to make this estimate small enough. 

We now proceed with our calculation, using w in place of w, i.e., —q in 
place of (3 2 . Partition sum (|33p in two pieces: 

Z(G;q,w) = £< + £>, 

where 

£< = £ h(B,-q)fi[ m and = £ -q) ^ . 

BCL:\B\<n+m BCL:\B\>n+m 

Set Q = (-l) n q 2n + m + 1 ( q - \) m ( q - 2) n , and note that Q ^ 0. We'll show: 
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1. If \B\ < n + m then h(B, -q) = 0. 

2. If \B\ = n + m therU 



h(B,-q) 



Q, if (V, T U B) is bridge connected; 
0, otherwise. 



3. The set {B : \B\ = n + m and (V, T U B) is bridge connected} is in 1- 
1 correspondence with the set of solutions to the instance of #3-D 
Matching. 



Observations HH3] entail 



n+m 



where ./V is the number of solutions to the #3-D Matching instance. On 
the other hand, E > is crudely bounded as follows: 



\2>\= £ KB,- q ) q \ B \(^\ 

BCL:\B\>n+m \ * / 

< J2\h(B,-q)qW(p 1 /q) n+m+1 \ 
BCL 

< 2 \V\q\ T ^ V \q\ L \\p 1 /q\ n + m+1 

< (2q) 6n+7m+1 \/3 1 /q\ n+m+1 . 



Let Q = q n+m Q. Now, setting e (the bound on \@i/q\) so that 



(2qf n+7m+1 e<±\Q\, 



we have 

Z(G;q,w) 



iV 



iV 



> 



n+m 



Then, according to ([S3]) , setting 

(3n + 2m)(2g) 6n+7m+1 5 < §Q/^ +w 



Q(A/q) 



n+m 



1 

< -. 



(36) 



ensures 



Z(G;q,w) - Z(G;q,w) 



QPx 



n+m 



1 

< -. 



(37) 



^Bridge connected is a synonym for 2-edge-connected, i.e., connected and having no 
bridge, which is an edge whose removal would disconnect the graph. 



27 



Combining (j36l) and (|37|) yields the required estimate 



Z(G;q,w) 
Qff +m 




It remains to verify the three observations. Suppose the graph H = 
G\B/B contains a bridge e. Then 



where the — q/q comes from including e in A, which gives a weight of —q but 
reduces the number of components by one and the 1 comes from excluding e 
from A. The tree (V, T) has 2(n + m) leaves, so if \B\ < n + m there are at 
least two vertices in (V, TUB) of degree one. The unique edge e incident at 
either of these vertices is a bridge, and is a member of T; it is clearly also a 
bridge in H. This deals with Observation [TJ 

Suppose B C L is a set of link-edges of size n + m such that (V, TUB) 
is bridge connected. Every leaf of (V, T) must have some edge of B incident 
at it, and hence exactly one. Call such a B a pairing. If i? is a pairing then, 
for every triple (wi,Xj,y k ) £ M, either (i) {xj, a^}, {a ijfe , G 5 and 
{a ijk ,a ijk } £ B, or (ii) {a ijk , a ijk } £ -B and {xj, a ijk }, {a ijk , y k } (£ B. Let 
M' be the set of triples of type (i). By counting, \M'\ = n. So there is a 1-1 
correspondence between pairings B and sets M' C M containing n triples. 
We will now argue that, under this correspondence, bridge-connected graphs 
(V, TUB) are associated with solutions to #3-D Matching and vice versa. 

On the one hand, if M 1 covers all of W U X U Y } then it is easy to 
check that every edge in (V, TUB) is contained in a simple cycle of the 
form (r,Xj,aij k ,Wi,r) or (r,y k ,a ijk ,Wi,r) for some triple (wi,Xj,y k ) £ M', 
or a cycle of the form (ittj, a^k, a^k, Wi) for some triple {wi,Xj,y k ) £ M\M', 
Conversely, if (V, TUB) is bridge connected then in particular B is a pairing, 
which immediately implies that every element of X and Y is covered by some 
triple in M' . But also every Wi must be covered, since the only way to avoid 
{r, Wi} being a bridge is to have either {aij k , Xj} £ B or {a,ij k ,y k } £ B, for 
some j, k (and hence, in fact, both). This is Observation [3j 

Finally to Observation [2j If (V, TUB) is not bridge connected then it 
has a bridge e which is necessarily a tree edge. (The link edges join leaves 
of the tree (V, T), and hence cannot be bridges.) We have already seen 
that the existence of a bridge implies h(B, —q) = 0. So suppose (V,T U B) 
is bridge connected, and let M' be the corresponding 3-d matching. Then 
graph H = G \ B/B may be described as follows. 

For each triple t = (wi,yj,z k ) £ M', denote by H t = (Vt,Et) the graph 
with vertex set 



Z(H;q,w) 



A'CE{H)-{e} 




w 




V t = {r,w i: Xj,y k } U { 



dij'k' : (iVi,Xj>,y k >) £ M \ M'} 
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and edge set 



E t = {{r.Wi},^,^},^,^},^,^},^,^}} 
U {{wi,a ij/k/ } : (wi,Xji,y k >) G M\ M'}. 

The edges with endpoints of the form aijiy have multiplicity two, the others 
multiplicity one. Observe that VtH Vf = {r} for distinct triples t' ^ t. (This 
is a consequence of M' being a 3-d matching.) The graph H is obtained 
by taking the union of all Ht and identifying vertex r, so Z(H ; q, w) = 
Q l n Tit Z(Hf, Q,w). Each of the multiplicity-two edges (m — n of them) 
contributes a factor q{q — 1), which is non-zero by assumption. That leaves 
us with n copies of K4 minus an edge. Each of those contributes a factor 
— q 4 (q — l)(q — 2), which again is non-zero, by assumption. Putting it all 
together, 

Z(H; q, w) = q X - n [-q\q - l)(q - 2)] n [q(q - l)] m - n 
= {-l) n q 2n+m+1 {q - l) m (q - 2) n . 

Finally, we need to relate our conveniently chosen weights, j3\ and P2, 
to the actual ones, a\ and «2- This is done as in the proof of Lemma [T3l 
In particular, we choose k\ and &2 satisfying (q/a\ + l) kl — 1 > 1/e and 
\(q/a 2 + if 2 1 < 6. Now G' is formed as in the proof of Lemma [131 except that 
/c-stretches are used in place of fc-thickenings (according to Equation I14p . 
As before, k\ = 0{m) and k 2 = 0(m 2 ), so the construction is polynomially 
bounded. □ 

Using Lemma we can now prove Theorem [3l 
Theorem [3] Suppose (x,y) £ Q 2 satisfies q = (x - l)(y - 1) ^ {0,1,2}. 
Suppose also that it is possible to shift the point (x, y) to the point (x 1 , y') 
with x' £ [-1, 1], and to (x" , y") with x" G (-1, 1). Then there is no FPRAS 
for Tutte(x, y) unless RP = NP. 

Proof. The proof is analogous to the proof of Theorem [2l Note that none 
of y, y' and y" are equal to 1 since q 7^ 0. Let a = y — 1, a\ = y' — 1 and 
«2 = y" — 1- The constraints on x' and x" , together with {x 1 — l)(y' — 1) = q 
and (x" - l)(y" - 1) = q imply that x' - 1 = q/a x (£ [-2, 0] and x" - 1 = 
q/a2 G (—2,0). The proof is now exactly the same as the proof of Theorem [2] 
except that Lemma [15] is used in place of Lemma [13j □ 

4.5 The hyperbola H 2 in the halfplane y < — 1 

The following lemma will be used in the proof of Lemma \7\ 

Lemma 16. Suppose a\,a 2 G Q — {0} satisfy 2/a\ ^ [—2,0] and 2/a 2 G 
(-2,0). Then ^Perfect Matchings <ap MultiTutte(2; ai, a 2 ) 
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Proof sketch. The construction and proof of correctness are simplified ver- 
sions of those from the proof of Lemma [T5l so we provide only a sketch 
here. 

Suppose G = (V,E) is an instance of ^Perfect Matchings. For 
convenience, set n = \V\/2. Let G = (V,E) be the graph with vertex set 
V = VU{t} and edge set E = E U T, where T = {{t, v} : v G V}. Define 
w : E — > { /5i , /?2 } by w(e) = f3i if e £ E and 10(e) = ft in e £ T. As before, 
is small in absolute value, and ft is close to — q = —2; specifically, 
\Pi/q\ <eand |l + 2/ft| < 6. 

Following the now familiar path, 



Z(G;2,w)= ^^)2 K 



^ ^ W (BUC)2 K ( BUC ) 

5CBCCT 



where 



= 2 h{B,fh)pr, 

BCE 

h{B;(3 2 ) = M C| 2 K(BUC) (38) 

CCT 



Set Q = q n+1 (q — l) n = 2 n+1 . We will show the following observations. 

1. If \B\ < n then h(B, -2) = 0. 

2. If \B\ = n then 



h(B,-2) 



Q, if (V, TUB) is bridge connected; 
0, otherwise. 



3. The set {B : \B\ = n and (7, TUB) is bridge connected} is in 1-1 
correspondence with the set of solutions to the instance of ^Perfect 
Matchings. Specifically, (V,T U B) is bridge connected iff B is a 
perfect matching in G. 

Thus, for e,S sufficiently small, |2-( n+1 )/3f n Z(G; q, w) - N\ < \, where N 
is the number of perfect matchings in G. The proof is completed exactly as 
before. 

It remains to justify the three observations. For Observation 1, note 
that if \B\ < n then (V,B) contains an isolated vertex. Consider the factor 
contributed to h(B, —2) from the edge connecting this vertex to t. The 
contribution is —q (for including the edge) plus q (for excluding it, and 
hence adding a component), which is 0. Observation 3 is self-evident. Using 
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Observation 3, we can establish Observation 2 as follows. Suppose that B 
is a perfect matching. Then 

h(B, -2) = 2((-2) 2 + 2(-2) + 2)" = 2 n+1 , 

where the first 2 comes from the component containing t, and for each of the 
n edges in the matching, the (— 2) 2 comes from including both edges to t, 
the 2(— 2) comes from the two ways to add one of the edges to t, and the 2 
comes from excluding both edges to t, which adds a component. 

□ 

We can now prove Lemma [7J 
Lemma Suppose (x, y) is a point with {x — l)(y — 1) = 2 and y < — 1. 
Then Tutte(x, y) =ap #Perfect Matchings. 

Proof. To show #Perfect Matchings <ap Tutte(x, y) use thickenings 
as in the proof of Corollary [5] to shift (x,y) to a point (x',y') with x' £ 
[—1,1] and to a point (x",y") with x" G (—1,1). Then follow the proof 
of Theorem [3] to reduce MultiTutte(2; ai, a 2 ) to TuTTE(:r,y). Finally, 
Lemma reduces ^Perfect Matchings to MultiTutte(2; a\, 02). 

We now show Tutte(x, y) <ap #Perfect Matchings. Using the 
Fortuin-Kasteleyn representation of the Potts model, Z(G; 2, y — 1) is equal 
to the partition function of the Ising model in which every edge has weight 
y — 1. That is, 

Z(G;2,y-l) = £ y™^, (39) 
«V(G)-*{-l,l} 

where mono(ff) denotes the number of monochromatic edges in the map- 
ping a. (See (T7J (2.7), (2.9)] for a justification of this identity.) We will 
assume without loss of generality that the graph G has no loops. It is clear 
from (f39j) that a loop merely introduces a factor of y. 

Let n = \V{G)\ and m = \E(G)\. Let v = (y - l)/(y + 1). Now let G' 
be the graph constructed from G by replacing each vertex of degree i > 4 
as follows. Suppose that the neighbours of vertex v in G are w\, . . . ,wi. 
Then replace v with a path of £ — 2 new degree-3 vertices i^, . . . 
The edges (^2, ^3), (^3, V4), . . . , (1^-2, V£-i) will be called "supplementary" 
edges of G' . The edges (wi,v), (u>2,v), . . . , (u;£_i,u), (w£,v) correspond to 
edges (wi,V2),(w2,V2),(w3,vz),...,(wt-. 1 ,V£-i),(we,Vi- 1 ) of G' . We will 
call these edges "primary", because they correspond to the original edges 
of G. 

Then let G" be the graph constructed from G' by replacing each vertex of 
degree 2 as follows. Suppose that the neighbours of vertex v in G' are w\ and 
W2- Then replace v with two new vertices v\ and V2 and replace the edges 
(wi,v) and (^2,^) with the path (wi, v±), (v 1, f 2 ), (w 2 , W2) in which (wi,f 2 ) 
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is a supplementary edge of G" and the edges {w\,v\) and (t>2,u>2) are "pri- 
mary" edges of G" . Finally, if v is a degree-3 vertex of G' with neighbours 
wi, W2 and W3, replace v with the three vertices vi,V2,vs. Add supplemen- 
tary edges (vi,V2), (v2,vs), (v3,vi). Replace the edges (wi,v), (w2,v) and 
(w3,v) of G' with edge (w\,vi), (102,^2) and (^3,^3) of G" , making edge 
(wi,V{) primary if and only if (vii,v) was primary in G' . (Once again, the 
primary edges of G" correspond to the original edges of G.) 
Fisher has shown [7J (10)] that 

(\ m -. 

where the sum is over perfect matchings X of G" and the product is over 
primary edges e of G" that are in the perfect matching X. 

Now let n\ and ri2 be positive integers such that 1/v = n\jn%. Let H 
be a graph consisting of n\ parallel edges from a vertex u to a vertex a and 
n,2 parallel edges from the vertex a to a vertex b and a single edge from b to 
a vertex v. Let KA be the set of matchings of H which match both a and b. 
There are n\ matchings in A4 in which a is matched with u. All of these 
match both u and v. There are 712 matchings in A4 in which a is matched 
with b. These do not match u or v. There are no other matchings in Ad. 

Construct G from G" by replacing every primary edge (u, v) of G" with 
a copy of H. Then the expression Ylx lie \ * n Equation (fiOj) is equal to the 
number of perfect matchings of G divided by n™. So if we could approximate 
the number of perfect matchings of G, we could approximate Z(G\ 2, y — 1). 

□ 

Remark. The construction in the reduction from the problem Tutte(x, y) 
to the problem ^Perfect Matchings relies on the fact that y — 1 and 
y + 1 have the same sign (so n\ and ri2 are both positive integers). The 
same reduction would apply for q = 2 and y > 1 but this is ferromagnetic 
Ising, and we already have an FPRAS, due to Jerrum and Sinclair |12j . 

5 #P-hardness 

In Section 12.3} we noted that if x > and y > then Tutte(x, y) is 
in #Pq>, so there is a randomised approximation scheme for Tutte(x,?/) 
that runs in polynomial time using an oracle for an NP predicate. Here 
we show that it is unlikely that Tutte(x, y) is in #Pq for all x and y. In 
particular, we identify a region of points (x, y) where y is negative for which 
even approximating Tutte(x, y) is as hard as #P. Specifically, we prove the 
following. 

Theorem 17. Suppose (x,y) is a point withy E (—1,0) and (x — l)(y — l) = 
4. Then there is no FPRAS for TuTTE(x,y) unless RP = #P. 
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5.1 The Potts model 

For a positive integer q and a y 6 Q, and a graph G = (V, E), let 
P(G;q,y)= y mono{(7) , 

a:V->{l,...,q} 

where mono(cr) is the number of edges in E that are monochromatic under 
the map a. P(G; q, y) is the partition function of the q-state Potts model 
at an appropriate temperature (depending on y). The region y > 1 is 
"ferromagnetic" since like spins are favoured along an edge, the region < 
y < 1 is "antiferromagnetic" , and the region y < is "unphysical" [T7j. It 
is known that the g-state Potts model coincides with the Tutte polynomial 
when q is a positive integer. In particular (see (fTU|) and [TTJ (2.9)]), 

T(G; x, y) = (y- l)" n (x - 1)" k(£;) P(G; q, y), 

where q = (x — l)(y — 1). 

In the rest of this section, we suppose that we have an FPRAS for 
P(G; 4, y) for a point y E (—1,0) and we show how to use the FPRAS to 
solve a #P-hard problem (counting proper 3-colourings of a simple graph). 

First, we establish some notation. If G is a graph with designated vertices 
a and b and a and /9 are values in {1, . . . , q}, let P(G; q, y \ a(ab) = af3) 
denote the contribution to P(G; q, y) due to colourings a with a(a) = a and 
a(b) = (3. 



5.2 The building blocks 

Fix y £ (—1,0). Suppose that n is the number of vertices of a graph G. Let 

2 

M be a rational number in the range 1 < M < 3 n and let e = 2~ n . In 
this section, we will show how to construct a graph Hm with two designated 
vertices, a and b, so that 

-1 P(ff M ;4,y |<7(afr) = 11) -1 

M - P(H M -A,y | a(a6) = 12) ~ M ' 1 J 

As a building block, let Pf be an £-edge path. Let fi denote P(Pf,A, y \ 
a(ab) = 11) and let ai denote P(P^;4, y \ a(ab) = 12). These satisfy the 
recurrences ft = yfe^i + 3a^_! and = f^-i + (2 + y)a^_i with fi = y and 
ai = 1. The solution to these recurrences, for I > 1, is given by 

and 

a e = \(3 + y) i -\(y-l) e . 
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Thus, 

at _ (3+ y y-(y-iy 4(y-iy 

h (3 + y y + 3(y-iy (3 + y y + 3(y-iy 

Recall y > -1 and let 7 = ((3 + y)/{y - I)) 2 > 1. For every positive integer 
j, let 5j = 1 — a,2j/f2j- Note that 

< 7 ~ j < 5j < 4 7 - J '. (42) 

Also, fzj/aQj = 1/(1 - 

Given y, choose a positive odd integer k so that 

\y\ k < ± < \y\ k - 2 . (43) 

Now, let t be the smallest integer such that St < eM. For j € {1, . .. ,£}, 
choose a natural number Ay so that 

i-± 1 1 l jl l 1 

|y ' fc 2 (1 - W r (1 - ^ ~ M < ly|fe n (1 _ 5 r )k r (i _ ^jvPi • ( 44 ) 

Now i?M is formed by joining a number of paths, all with endpoints a 
and 6. To form Hm, take A; paths of length 1 (i.e., edges). Also, for every 
j G {1, . . . , t}, take Ay paths of length 2j. So 

P(H M ;A,y\a(ab) = 11) _ ' 1 

P(H M ; 4, y | a{ab) = 12) " 12/1 ^ (1 - 5 r )^ ' 

and this is at least — 1/M and at most — (1/M)(1 — St), which implies Equa- 
tion flu}. 

Now we consider the size of Hm- Equation (|43p implies that k = 
O(logM) = 0(n). Also, t = 0{n 2 ) by (02}. 

How big can Ay be? By (gU) we have {1 _] )k] < , so (1 - ^)^' > 

1 — (5j_i . (5j is decreasing in j, so without loss of generality, we'll deal with 
those j such that Sj < 0.7 (the values of kj corresponding to smaller values 
of j are just constants). Then 

(1 - Sjf 1 ** < 0.15 < (1 - Sj-i) 1 /**- 1 , 

so 

(1 - 63)™*-"-'** < (1 - Sj-xfi- 115 '- 1 = 1 - 6j-x, 
and therefore kj < 26j-\/5j = O(l). 
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5.3 The construction 



We use the notation from Section [S~2l Let r be the smallest even integer such 
that \y\ r < eAT n . Construct G' from the simple graph G (the graph we wish 
to 3-colour) by replacing every edge (u, v) of G with a bundle of r parallel 
edges with endpoints u and v. (That is, we perform an r-thickening on all 
edges.) Add two new vertices, a and b. Connect a to every vertex in G by 
a bundle of r parallel edges. Similarly, connect b to every vertex in G. 

Let n denote the number of vertices of G and m denote the number of 
edges of G. Recall that P(G; 3, 0) is the number of proper 3-colourings of G. 
Then, 

P(G; 3, 0) < P(G'; 4, y \ a(ab) = 11) < P(G; 3, 0) + 4 n |y| r , 

so 

P{G\ 3, 0) < P(G'; 4, y \ a{ab) = 11) < P(G; 3, 0) + e. (45) 

Similarly, 

P(G; 2, 0) < P(G'; 4, y \ a{ab) = 12) < P{G; 2, 0) + e. (46) 

Let Gm be the graph constructed from G' and Hm by identifying vertex 
a in G' with vertex a in Hm and similarly identifying vertex b in G' with 
vertex b in Hm- Let 

P{G M -A,y) 



P(H M ;4,y | a{ab) = 12) 

and note (from the previous section) that the quantity P(Hm', 4, y j a(ab) = 
12) in the denominator is positive. Now 

P(G M -A,y | a(ab) = 11) = P(H M ;4, y \ a(ab) = ll)P(G';4,y | a(ab) = 11), 
and 

P(G M -A,y) =4P(G M ;4,y I <r(a6) = 11) + 12P(G M ; 4, y \ a(ab) = 12), 

so 

P(H M -A,y | <r(a&) = 11) p ,„, , , . . 
P{H M \%y | o-(ao) = 12) 
+ 12P(G';4,y | a{ab) = 12). 

Let £ = 5e3 n = o(l). By Equations flU}, gSJ, and gBJ, 

= "4 P( ^ 3,0) + 12P(G; 2, 0) + g M , 

where \qm\ < 
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We will restrict attention to graphs G which are bipartite with at least 
4 vertices. Note that it is #P-hard to count the proper 3-colourings of a 
bipartite graph. For example, [3J Section 6] observes that this is the same 
as counting homomorphisms from a general graph to the cycle Cq, which is 
shown to be #P-hard by Dyer and Greenhill [6]. Also, for such a graph G, 
P{G; 2, 0) > and P{G; 3, 0) > 4P(G; 2, 0). 

Now, suppose that we had an FPRAS for approximating P(Gm', 4, y)- A 
call to the FPRAS gives us the sign of Ym ■ 

Let G be a bipartite graph with n > 4 vertices. Let zi = 3~ n and z u = 1. 
Then we have an interval [zg, z u ] with Y\i z > and Y\i Zu < 0. Use binary 
search to bisect the interval until we have an interval [zg, z u ] with Y\i Zl > 0, 
Yi/ Zu < 0, and z u — zi < e. (This takes at most re 2 bisections since, after 
j bisections, z u — z^ < 2~- J .) 

Since Y\/ Zl > 0, we have 

^ 3P(G;2,0) f 



Similarly, since Y\i Zu < 0, 
z u > 

So 



P(G;3,0) 4P(G;3,0)" 

3P(G;2,0) g 
P(G;3,0) 4P(G;3,0) 



3P(G;2,0) g 3P(G;2,0) g 

P(G; 3, 0) 4P(G; 3,0) - z « - ^ + e - jP ( G; 3; ) + 4P(G; 3, 0) £ ' 

Now the point is that only one real number in the specified interval for 
z u is of the form 3ni/n2 where n\ is an integer in {1, . . . , 2 n } and is an 
integer in {1, . . . , 3 n } (since e and £ are so small) so the value of z u allows 
us to compute 3P(G; 2, 0)/P(G; 3, 0) exactly, and since P(G;2,0) can be 
computed exactly, this gives us P{G; 3, 0), thus counting proper 3-colourings 
of G. 
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